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ABSTRACT 


The  equations  and  a  Fortran  program  to  calculate  supersonic  and 
hypersonic  aerodynamic  heat  transfer  rates  and  transient  temperature 
distributions  for  spherical  leading  edges  and  flat  plate  surfaces  are 
presented  in  this  report.  The  missile  skin  is  composed  of  one  to  three 
different  slab  materials  and/or  thin  wadi  combinations  for  flight  trajec¬ 
tories  or  wind  tunnel  conditions.  The  Fortran*  program  is  written  for 
the  IBM  1620  40K  digital  computer. 


CONTENTS 


ABSTRACT, 


1.  INTRODUCTION 


ANALYSIS 


a.  Stagnation  Regions . . 

b.  Flat  Plate  Regions . . 

c.  Time  Increment  and  Material  Properties  . 

d.  One  Dimensional  Temperature  Distribution 

e.  Specific  Heat  Ratios  for  Air  .  . . 

f:  Flight  Environment  ........... 


3.  CONCLUSIONS  . 
LITERATURE  CITED 


Page 

ii 

1 


1 

10 

12 

13 

17 

18 


18 


19 


Appendix  *A 

FORTRAN  PROGRAM  AND  IT5i  USAGE 

1.  Fortran  Program  Statements . 21 

2.  Input  Format  . .  27 

3.  Input  Comments .  28 

4.  IBM  1620  Operating  Instructions .  29 

5.  Example  Runs  for  Sphere,  Flat  Plate,  and  Cone  ...  31 

6.  Input  Data  for  Examples  . .  31 

7.  Output  Data  for  Examples . 34 

Appendix  B.  1959  ATMOSPHERIC  PROPERTIES .  39 


SYMBOLS 


41 


ILLUSTRATIONS 


r.t  ' 

K. 


Important  Variables  Affecting  the  Aerodynamic  Heat 
Transfer  Coefficient  for  a  Spherically  Blunted  Leading 

Edge  Surface . 

Stagnation  Point  Velocity  Gradient  .  . . 

Laminar  Loading  Edge  Skin  Friction  Proportionality 

and  Velocity  Gradient.  . . 

Turbulent  Leading  Edge  Skin  Friction  Proportionality 

and  Velocity  Gradient . 

Laminar  Heat  Transfer  Distribution  , . 

Aerodynamic  Heat  Transfer  Variables  for  Flat 

Plates  or  Cones . 

Multi -Slab  Materials . . 

Thin  Wall  Followed  by  Multi-Slab  Materials . 

Multi -Slab  Materials  Followed  by  a  Thin  Wall  .  .  .  .  , 


(M 

:,"(V 

IT  . 


iv£$ 


Table 


Page 


I  ARDC  1959  Atmospheric  Properties 


Figure 


iv 


1.  Introduction 


A  general  purpose  Transient  Temperature  Aerodynamic  Heat 
Transfer  IBM  1620  Digital  Computer  Program  for  supersonic  and  hyper¬ 
sonic  flight  speeds  is  described  herein.  This  computer  program  con¬ 
siders  spherically  blunted  leading  edges  and/or  flat  plate  surfaces.  One 
dimensional  temperature  distributions  through  a  missile  skin  composed 
of  one  to  three  different  slab  materials  or  a  thin  wall  material  followed 
or  preceded  by  one  or  two  different  slab  materials  is  available.  The 
required  flight  environment  is  either  a  trajectory  input  based  on  the 
ARDC  1959  atmosphere  or  constant  altitude  and  local  flow  properties 
(wind  tunnel  conditions). 

2.  Analysis 

a.  Stagnation  Regions 

Aerodynamic  heat  transfer  coefficients  for  spherically 
blunted  leading  edge  surfaces  are  separated  into  two  regions.  For  non- 
dissociated  gas  properties,  corresponding  to  flight  speeds  up  to 
6000  ft/ sec,  the  external  aerodynamic  heat  transfer  coefficients  for 
laminar  and  turbulent  boundary  layers  are  developed  in  this  report.  For 
dissociated  gas  properties,  an  approximation  to  the  exact  stagnation 
point  heat  transfer  rate  soluuon  of  Fay  and  Riddell1  and  Detra  and 
Hidalgo2  is  used. 

( 1)  Nondissociated  Aerodynamic  Heat  Transfer  Coefficient 
The  important  variables  affecting  the  ae  rodynamic  heat  transfer  coeffi¬ 
cient  for  a  spherically  blunted  leading  edge1  surface  are  shown  in 
Figure  1. 


Figure  1.  Important  Variables  Affecting  the  Aerodynamic  Heat  Transfei 
Coefficient  for  a  Spherically  Blunted  Leading  Edge  Surface 


(a)  Laminar  Boundary  Layer.  A  modified  Reynolds 
analogy  for  flow  with  constant  thermal  and  transport  properties  through 
the  boundary  layer  for  spherical  and  cylindrical  surfaces  is  used. 
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Nuinc  Cf  "0.6 
Res  Prs  =  2  Pr® 


(1) 


•^^-=<U-ResPr,»-4  .  (2) 

Measurements  of  local  skin  friction  coefficients  on  spherical  and 
cylindrical  surfaces’  indicate  the  normal  laminar  flow  correlation  for 
zero  pressure  gradient  flow  may  be  applied  provided  the  constant  of 
proportionality,  using  local  Reynold's  number, is  considered  to  vary 
with  location  along  the  surface: 


Cf  = 


-s fTCeZ 


(3) 


The  factor  fi  of  Equation  (3)  varies  from  1 . 526  at  the  stagnation  point 
to  approximately  0.  664  at  a  position  90  degrees  from  the  stagnation 
point  for  the  sphere  and  from  1.  14  to  0.  664  for  the  cylinder.  Equation 
(2)  can  be  expanded  to 


H 


inc 


Y  =- 


[~P  a  Y  g  Y  ]  °'! 


0 


Pr, 


(4) 


Since  the  stagnation  point  value  of  Y  and  V8  are  zero  and  large  errors 
of  heat  transfer  coefficient  result  from  Equation  (4)  in  the  areas  close 
to  the  stagnation  point,  it  becomes  convenient  and  more  accurate  to 
define  a  term  p  as  given  by  Equation  (5): 


V3  =  PY  .  (5) 

The  value  of  p  varies  along  the  surface,  with  the  sphere  diameter,  and 
with  the  free  stream  Mach  number  and  temperature.  At  the  stagnation 
point,  application  of  elementary  calculus  yields  the  nondimen sional 

(dV\ 

—y  1  .  This  nondimensional  velocity  gradient  is  a 
function  of  free  strearivMacn  number  only.  A  ratio  of  P/Po  was  found 
to  depend  only  on  the  location  on  the  spherical  surface. 
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Substituting  pY  for  V8  and  ~-=—  for  pg  in  Equation  (4)  yields 

«-gl8 


Hinc*Y  = 


1  k  r  p8pY2T,sPi  o.4 

2  8  |_RgT » ^sj  8 


By  dividing  Equation  (6)  by  Y,  multiplying  by  n/D,  then  multiplying  tt>.e 
right  hand  side  by  ^p0/p0  and  's/Voo/Voo,  the  equation  may  be  written: 


where 


•  Vd  =  0.  5  2,ZtZ>a  . 

«.g 


M  /  r-  T  f"  12»5i 

Z,=^H[14  +  m^]  [°-,s»(T-iS?-)]  j  .  «8> 


Z,  =  kPr0’4 


Equation  (7)  has  been  developed  for  constant  thermal  and  transport 
properties  with  Z3  evaluated  at  local  conditions.  For  an  appreciable 
variation  of  temperature  within  the  boundary  layer,  a  reference  temper¬ 
ature3  (T#)  has  proven  to  give  excellent  aerodynamic  heat  transfer 
coefficients: 


T*  =  Ts  0.50  +  0.  039  Mfl 


']  +0 


.50  T„ 
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The  properties  of  Z5  in  Equation  (10)  required  to  be  evaluated  at  T*, 
are  indicated  by  the  asterisk  superscript  (*),  and  are  defined  by  the 
following  equations: 


When  T#<  1000eK 

,  0.  23791763  X  10'6  T* 

k*  - - 

T*  +  198.  6 

When  T*  2  1000°R 


k*  =  1 1 . 997  (i*, 


(12a) 


(12b) 


Where 


0.  249  X  10-*  T*  0,63 
~  8a 


(13) 


Pr*  =  P*  c9  Sa/k’ 
and  Cp*  =  f  (T*),  as  defined  by  Equation  (56). 
Then 


(14) 


Hn/5  =  0.  5 


rv»p.i 
L  Rg  j 


0.5 


z.z2z,  . 


(15) 


Figures  2  and  3  present  Z(  and  . 

At  the  stagnation  point,  the  rji  -ren-e  temperature  becomes 

T0*  =  0.  5  JttqT  +  Tw]  for 


Mi 


-0.  0  , 


(16) 


where  the  total  temperature,  Tq-QT*  at  the  stagnation  point  is  the  otal 
temperature  of  the  free  stream. 


(b)  Turbulent  Boundary  Layer.  It  is  possible  to  have 
turbulent  boundary  layer  flow  over  some  portion  of  the  leading  edge 
and  a  method  for  the  aerodynamic  heat  frartsfer  coefficient  is  presented. 
The  basic  development  is  similar  to  the  laminar  boundary  layer  heat 
transfer  development. 
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The  skin  friction  coefficient  for  turbulent  boundary  layers  on  the  leading 
edge  can  be  expressed  as 


Cf  = 


0.2 


(18) 


The  leading  edge  geometry  has  only  a  minor  effect  on  the  value  of  the 
proportionality  constant  f2  as  compared  to  the  laminar  boundary  layer. 
The  final  equation  for  turbulent  leading  edge  aerodynamic  heat  transfer 
coefficient  is 


(19) 


(20) 


and 


Z,= 


1/3 

(TV5)  0,8 


k*Pr* 


(21) 


(22) 


The  e?cponents  of  Equation  (14)  reflect  the  basic  changes  in  the  skin 
friction  correlation.  Figures  2  and  4  show  variations  of  Zj  and  Z2  . 

The  term  (Y/D)0*6  will  cause  a  maximum  aerodynamic  heat  transfer 
coefficient  away  from  the  stagnation  point  on  a  given  surface  and 
theoretically  a  value  of  zero  at  the  stagnation  point.  In  reality,  the 
stagnacion  point  flow  is  laminar  and  thus  the  aerodynamic  heat  transfer 
coefficient  will  not  become  zero.  This  turbulent  analysis  is  not  included 
in  the  Fortran  program. 


(c)  Approximate  Pressure  Distribution.  The  reference 
temperature  ,  T$,  requires  the  local  Mach  number  and  local  temper¬ 
ature.  The  basic  aerodynamic  heat  transfer  equations.  Equation  (7) 
for  laminar  boundary  layer  and  Equation  (19)  for  turbulent  boundary 
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layers  ,  require  local  pressure  value.  A  modified  Newtonian -Prandtl 
Meyer  pressure  ratio,  Ps,  as  a  function  of  angular  position,  9  is  used. 


where 


P  -  p’  ps  pTOT 

~  3  "  p  — ^ 

i  O  i  ^oo 


=  1  -  0.  957  sin2  0  , 

PTOT 


The  local  Mach  number  is  determined  from  the  following  equations: 


Mife)' 


PTOT  -  1 1 (“Y  +  1) 


_ y_ 

2  y  hloo 


+ 1  1W-1 

‘-(y-DJ 


The  local  temperature  Ts  is  obtained  from: 


*  1 


a 


(y  -  1)  Moo2* 
(Y  -  1)  Ms2_ 


for  surface  positions  away  from  the  stagnation  point.  At  the  stagnation 
point,  the  reference  temperature,  T*,  of  Equation  (11)  does  not  require 
local  temperature. 

(2)  Dissociated  Air  Aerodynamic  Heat  Transfer  Rates. 

An  approximate  equation  for  the  exact  stagnation  point  aerodynamic 
heat  transfer  rate  for  flight  velocities  greater  than  6000  ft/sec  is 
presented.  2 

QwWi^  =  865(^\|TH -  fe-0--.--?:  .  (28) 

*  \10  /  If*sea  level  ho  _  hw,00 


where 


hQ  =  stagnation  enthalpy  =  6006  Tw  +  0.  5  Vo,2  (28a) 

hw=  enthalpy  at  Tw,  "R  =  778  ga  cp  Tw  (28b) 


and 


f 


£ 

£ 

f 


h„,  =  enthalpy  at  300,  °k  =  3244100. 

"  300  ' 


(28c) 


For  variation  of  laminar  heat  transfer  rates  around  the  spherical 
blunted  leading  edge,  the  Lee's  ratio  of  heat  transfer  rate  to  stagnation 
point  heat  transfer  rate3  QR  is  presented  in  Figure  5. 

Turbulent  leading-edge  boundary  layer  heat-transfer  rate  analyses 
are  not  included  for  the  hypersonic  flight  speeds. 

b.  Flat  Plate  Regions 

The  aerodynamic  heat-transfer  coefficients  for  laminar 
and  turbulent  boundary  layers  over  a  flat  plate  and/or  cone  surface  were 
developed  in  detail.4  Aerodynamic  heat  transfer  variables  are  illustrated 
in  Figure  6 . 


Figure  6.  Aerodynamic  Heat  Transfer  Variables  for 
Flat  Plates  or  Cones 
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T  = 


rtot/NLAY 


(37) 


gf- 

W 

sc 

IT. 

& 


*5 


& 
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At  = 


0.  5pCpT2 

k  +  VjT 


(38) 


where  V  =  10  for  first  material,  estimate  for  maximum  value,  and 
V  =  0  for  following  materials.  The  time  increment  should  be  approx¬ 
imately  the  same  for  all  the  materials  used. 


used. 


Equation  (38)  is  solved  for  each  material  and  the  smallest  value  is 


Other  required  material  functions  are: 


=  At/(p  C  t) 


(39) 


and 


(40) 


(41) 


where  m  is  the  number  of  the  material. 


d.  One  Dimensional  Temperature  Distribution 


below. 


The  basic  heat  balance  for  a  multi-material  skin  is  developed 


NL  = 


NMAT  = 
M  = 

■V  _ 

T’  = 


Total  number  of  layers  for  all  materials  plus  end 
point  (limited  to  15  in  program). 

Number  of  local  point,  from  1  to  NL. 

Total  number  cf  materials  (limited  to  3  in  program). 
Number  of  material,  from  1  to  NMAT. 

Temperature  for  each  point  at  the  present  time  step. 
Temperature  for  the  point  at  the  previous  time  step. 


13 


The  temperature  increment,  AT,  for  any  point  is  the  temperature 
difference  between  the  present  and  previous  time  steps. 

AT  =  T  -  T'  (42) 

Then  temperature  increments  between  local  layer  and  other  layers 
at  previous  time  step  are  defined. 

=  T'l+i  -  T'L  (43) 

D2*T«l_x  -  T't,  (44) 

D3  =  T'3  -  T'j  (45) 

For  all  cases,  heat  in  =  heat  out  +  heat  stored,  or 
Tin  =  Tout  +  ^stored 

(1)  Multi -Slab  Materials.  The  multi-slab  materials  are 
shown  in  Figure  7. 


interface  between  materials 


point 


Figure  7.  Multi-Siab  Materials 

(2)  Heat  Balance  at  the  Air  Flow  Side  of  Slab.  Heat 

in  =  qjn  =  H  (Trec  -  T})  =  Aerodynamic  heat  transfer  rate.  Heat  out 
is  the  heat  transfer  to  surrounding  atmosphere  plus  the  heat  transfer 
to  the  next  thickness. 


Tout  =  €<r  <Tr  +k  (TJ  -  Tz)'t 


^>1 


> -^S*"**/-*?* i&gZfuS  ’  ->J.w-vt*3,-afc*W^jS<r«SS5Z,^^RLfSS^&ac3fej 


Heat  stored  is  the  heat  remaining  in  the  first  thickness: 
^stored  =  0.  SAT^/Fj 


T  =  T 
rec  A  s 


[l  +  0.  5  (Vs  -  1)  RMS2  ]  . 


i ■  n 


For  a  Sphere:  Ys  =  1.4  . 

For  a  F.  P.  or  Cone:  Ys  is  defined  by  Equation  (57). 

Equation  (28)  defined  Q. for  a  sphere  when  the  local  velocity  is 
greater  than  6000  ft/sec,  but  for  other  cases, 

Qw  =H(Trec  -  Tj)  (4 

Q  -  Qw  -  e<r  (T  i  / 1 00)4  ,  (4 

With  the  heat  balance  qin  =  qout  =<lstored»  outside  skin  transient 
temperature  for  a  slab  is 


■.=  T',+  2[F,Q+B,D,]  • 


(3)  Heat  Balance  Around  Interior  Points. 


qin  M^m-i  "  Tm) /Tm 


qStored=-^  (T-TO  . 

The  temperature  at  each  small  layer  within  the  material  is 


TL  =  Tl/  +  Bmj^Dj  +  D2J 


where 


\.y 


L  =  2, 3  .  .  .  (NL-1). 


■% 


For  a  thin  wall,  the  temperature  is  assumed  to  be  constant  through  the 
entire  thickness,  thus 


t2  =  t1 


(7)  Multi -Slab  Materials  Followed  by  a  Thin  Wall.  The 
multi- slab  materials  followed  by  a  thin  wall  are  shown  in  Figure  9. 


interface  between 
slab  and  thin  wall 


►TNMAT-i- 


rNMAT 


NL-i 


tNL 


-m  =  NMAT  -  1- 


•m  =  NMAT 


Figure  9.  Multi-Slab  Materials  Followed  by  a  Thin  Wall 


For  a  thin  inner  wall,  when  NLAY,.,. .  .  _  =  1: 

NMAT 


TNL-i  = 


Vi  Bm 


where  m  =  innermost  material  -  NMAT, 


TNL  =  TNL-i 


since  the  temperature  is  assumed  constant  through  the  thin  material. 


«'S8S®*S*?-  -yg  - - . - 


ssp^sSg 

^PH 


^J®p 

s-jS-SiisS 


e.  Specific  Heat  Ratios  for  Air 


A  correlation  for  specific  heat  ratio  for  air  is 


Where  T 


-  s 

= 

0 

Cp  = 

0.24 

= 

800 

0.  24 

= 

1700 

= 

0.  27 

= 

3000 

0.  29 

- 

5000 

0.  31 

=;• 

9000 

= 

0.  32 

= 

11, 700 

= 

0.40 

= 

14,400 

Cp  = 

0.46 

^iocal  —  Cp/  [cp  -  Rg/J  J 


where 


T  „„  ft  -  lb 

J  =  778  — — - x  ga 

Btu  Ba 


Flignt  Environment 


The  IBM  1620  digital  computer  program  has  the  ARDC 
1959  atmosphere  subroutine  as  an  integral  part  of  the  transient  aero¬ 
dynamic  heat  transfer  calculation.  Appendix  B  describes  this  subroutine. 
In  addition,  a  constant  altitude  and/or  constant  local  flow  properties  flight 
environment,  such  as  wind  tunnel  testing,  is  included  in  the  computer 
routine  for  flat  plate  or  cone.  The  symbol,  NCFIT,  determines  whether 
trajectory  data  or  a  constant  value  for  altitude  is  used.  If  NCFIT  is 
given  0,  constant  altitude,  local  Mach  number,  pressure,  and  temper¬ 
ature  are  given. 


3.  Conclusions 

The  aendynamic  heat  transfer  and  transient  temperature 
distribution  computer  program  described  in  this  report  provides  an 
economical  preliminary  design  capability  for  heat  transfer  analysis. 
Comparison  of  transient  temperatures  with  PERSHING  Ballistic 
Missile  flight  test  data5  and  with  more  sophisticated  aerodynamic  heat 
transfer  digital  computer  programs  indicates  very  good  agreement. 
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Appendix  A 

FORTRAN  PROGRAM  AND  ITS  USAGE 


1.  Fortran  Program  Statements 


1000 

1001 
10G2 
1003 
100^ 

1005 

1024 

1025 
1034 
1054 

1006 

1007 

1008 
100? 
1010 
1011 
1C  1 2 

1013 

1014 

1015 
1C16 
10  17 
loia 
101? 


TEMPERATURE  VS  TIME  FOR  LAYFRS  OF  SPHFP.E*  FLAT  PLATE.  OR  CONE 
COMPILED  ON  IBM  1620  40K  BEGINNING  AT  08200 
SW  1  ON  —  FOR  CARD  01JTPU7  AT  ALL  TIME  STEPS 

SW  2  ON  —  FOR  CARD  INPUT  OF  -  DTI ME.  DTOUT.  TLO.  THI 

S W  2  OFF  -  FOR  TYPEWRITER  INPUT  OF  —  DTIME.  DTOUT.  TLO.  THI 

SW  3  ON  -  TO  END  RUN,  OUTPUT  LAST  STEP.  BRANCH  TO  NEXT  CASE 

TBODY  =  TYPE  OF  BODY 

=  1.0  FOR  SPHERE'  —  NCFIT  =■  ?. 

-  2.0  FOR  FLAT  PLATE  —  NCFIT  =5  OR  0 

'=  3.0  FOR  CONE  —  NCFIT  =  5  OR  0 

NCFIT  =  NO.  OF  CURVES  FITTED  — 

=  f>  FOR  F.P.  OR  CONE  --  GIVEN  CONSTANT  LOCAL  M,  P,  T.  ALT 
=  2  OR  5  —  CUSyES  FITTED  ARF  - 

1. ‘  ALTITUDE  VS  T I MF 

2.  VELOCITY  VS  TIME 

3.  LOCAL  M  VS  FREE  STREAM  MACU  NUMBER 

4.  PS/PF  VERSUS  FREE  STREAM  MACK  NUMBFR 

5.  TS/TF  VERSUS  FREE  STREAM  MACH  NUMBER/ 

-  TH=POS I T I  ON  ANG.  RNORY  =  NGSE  RADIUS.  QRRET=  Q' RATIO 

CONE  --  TH=IDEMTIFYING  ANG.  RNORY=SURFACE  L»  QRRET  =  RET 
NUMBER  OF  MATERIALS  —  1 .  2 .  OR  3 


SPHERE 
F.P.  OR 
NK  AT  = 
NL 


=  TOTAL  NUMBER  OF  SKIN  LAYERS  f  END  POiNT  — 
SUBSCRIPTS— S  FOR  LOCAL— F  FOR  FRFE  STREAM 


(1H 

(1H 

I1H 

(1H 

(1H 


FORMAT 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAi 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
DIMENSION 
DIMENSION 


) 

42H 
46H 
40H 
2  OH 


LAYERS 

LAYERS 

LAYERS 


TEMPERATURE  VS  TIME  FOR 
TEMPERATURE  VS  TIME  FOR 
TfMPFRATURE  VS  T IMF  FOR 
DTIME  FOR  EACH  MAT) 

I  4E15.8,  1:5) 

(1H  48H  SW1  OUTS  AL-L  STEPS— SW? 

(1H  18H  -  SW  3  ENDS  RUN) 

I1H  45H  INPUT— DTIME. DTOUT.  FLO, T^5 — ON 
(  1H  22H  APP  1620  MIN  REQUIRED) 
i  115) 

(  IE  1 5 .5 »  115) 

I  49H  DT I ME=  DTOUT  = 

THET  a=  y* 

THET  A  =  RN= 

TMELT=  OEFF= 

MATERIAL  —  DTIME.  M,  K»  RHO. 

-  FOR  EACH  OUTPUT  TIME  STEP 

TiME  ALT.KfT 

QW  QDOT 

MACH  NO.  PRESSURE 
MACH  NO.  PRESSURE 


OF 

OF 

OF 


MAXIMUM  =  15 


SPHERE) 

FLAT  PLATE) 
CONE) 


CARD  INPUTS  4  TIME  VALS) 


4  LINES  3Y  TY1 


THI  =  ) 


UH 
( 1M 
<1H 
(  1H 
(  1H 
CH 
I  1H 
(1H 
(1H 
(  1H 
t  1H 


TLO= 

RET=) 

0R  = ) 

FMI SSlV! YY=  ) 

CP.  TAUT.  NLAY ) 

VEL, FT/SEC  RE) 

H  TAUAB  -  HEAT  COEFFS) 

TEMPERATURE  --  LOCAL) 
TEMPERATURE  —  FREE  STREAM) 
TEMPFRATURP  AND  COUNT  FOR  EACH  SKIN  LAYFR) 

TIME.  ALT.  VEL  -  AT  WALL  MELTING  POINT) 

APi 11 .3) >  TrMP  t 15 ) rTEMPL! 15 ) .CONIC . 15) .NUMI5 ) 

RHO  I  3 )  .CP  (3  )  .TAUC  3)  .03)  »F(3)  ,F(3)  .V(5) 


37H 
36H 
44H 
46H 
35H 
39H 
49H 
•  ?H 
48H 
4  3H 
42H 


READ  A  TMOSPlir-f<  I C  PROPFRT  i  £s  FROM  AP  TABLE  OF  APPENDIX  B 


DO  1  K= 1 »  1 1 

READ  1005.  ALT,  TF,  PF 

AOfK.l)  =  A..7  »  3.2808333 

API*.?)  =  Tf  *  1.8 

AP1K.3)  t  pP  »  .020^854 

PRINT  10?4 

PRINT  10?S 


(B-I) 
(B  -2} 
(B~3) 
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INPM  FOR  EACH  CASE  —  BE61NS  AT  STATEMENT  2 
PUNCH  1000 
PUNCH  1000 

RFAO  1007.  TBODY .  NCFIT 

READ  100L,  TH,  RNORY.  ORRET 

REAO  1005.  TMELT.  OE.  E 

Y  =  RNORY 
RET  =  ORRET 
PRINT  1000 
PRINT  1000 

IF  ITBODY-2.0  >  301  .  302.  303 
301  PUNCH  1001 
PRINT  1001 
PUNCH  1000 
PUNCH  1010 
OR  =  ORRET 
SQRD  =  I2.*RN0RY>**.5 

22  =  1.522036+TH»J1.S36913F-03+TH*{-3.701602E-04+TH*-2.688E-06> )  (9) 
PSVPT  =  1.0  -  .957*51 N(TH/57»29578)**2  (24) 

60  TO  304 
30?  PUNCH  1002 
PRINT  1002 
PUNCH  1000 
PUNCH  1009 
HX  =  1.0 
60  TO  304 

303  PUNCH  1003 

PRINT  1003 
PUNCH  1000 
PUNCH  1009 

HX  =  3.0  **  .5  (32) 

304  PUNCH  1005.  TH.  RNORY,  ORRET 

PUNCH  1000 

PUNCH  1C11 

PUNCH  1005,  TMELT.  QE.  E 

PUNCH  1000 
PUNCH  101? 

PRINT  1000 
PRINT  1004 
REAO  1006,  NMAT 

NL  =  1 

C  MATERIALS  LOOP  —  COMPUTE  MINIMUM  OTIME  FOR  EACH  MATERIAL 
V2  =  10. 

00  4  M=i ,NMA7 

READ  1005.  C(M)»  RHO(M),  CP<M),  TAUT,  NUM(M) 

VI  =  NUM(M> 

TAUIM)  =  TAUT/VI  (37) 

NL=NL+NUM(M) 

N!1M(M+1)  =  0.0 

F(M)  =  RHO(M)*CP(M)*TAU(M» 

T  EMP  ( M  i  =  CCMJ/TASMM) 

OTIMP  =  ,5*F(M)/(TFMP(M1+V?)  (38) 

V?  =0.0 

PRINT  1007,  OTIME.  M 

PUNCH  1007,  DTIMF,  M 

4  PUNCH  1005.  C ( M 1 ♦  RHO(«l.  CPtM),  TAUT.  NUMIM) 

IF  { SFNSF  SWITCH  ?  )  400,  401 
400  REAO  1005,  071M-,  OTOUT,  TLO.  TH 1 

CO  TO  40? 


401  PRINT  1034 

ACCEPT  1005 ♦  DTIME 

ACCEPT  1005*  0T0UT 

ACCEPT  1005.  TLO 

ACCEPT  1005.  TH I 

40?  NT  =  f THI-TLO) /0T1ME  +  1. 

T  1 620  =  NT/10  #NL/4 
PRINT  1034 
PRINT  1005.  T 1620 
PUNCH  1000 
PUNCH  1008 

PUNCH  1005.  DTIME.  DTOUT .  TLO. 
00  6  K=1 ,NMAT 

E(M)  =  OT I Mc /F I M ) 

P1M1  =  F(M)*TFMP(M) 

IFCM-1 15.6.5 

5  F ( M  5  =  TEMP !M) /TEMPI M-l ! 

6  CONTINUE 
NL 1  =NUM I  1 ) 

NL2=NUMI?) 

C  READ  INITIAL  TEMPERATURES  AND  COUNT 

00  t  L=1 .NL 

?  READ  1007.  TEMPLIL).  LI 

G  *  1.4 
GM1  =  G-l. 

GPl  =  G+l « 

G 1  =  1./GM1 

G2  =  ?./GMl 
RG  -  1716. 

GA  =  A?, 174 

0A8  =  0.0 

TAU  *.0  =  0.0 

ABOUT=0.0 

PUNCH  1000 

IF  (NCFIT1  9,  8,  9 

C  CONSTANT  LOCAL  VALUES  GlVFN 


8 

RFAO  1005.  ALT, 

VFL 

READ  10^5 ,  SM. 
PUNCH  1016 

Pc  . 

TP 

PUNCH  1005.  SM, 

PS, 

TS 

9 

PUNCH  1000 

PUNCH  1013 

PUNCH  1014 

PUNCH  1015 

IF  (NCFIT5  10.  11. 

10 

10 

PUNCH  1016 

PUNCH  1017 

11 

PUNCH  1018 

PUNCH  1000 

TIME  =  TLO 

N  =  0 

C 

TIME  etFP  LOOP  — 

N=C  OUNTFR 

1  10 

N  =  N  +  1 

Tw  =  tpm°L( 1 ) 

LI  =  1 

X  =  TIMF 

IF  (NCFIT1  12.  4 C »  1? 
CURVE  FITS  DATA  LOOP 
DO  05  I =l «NCF!T 


(A  - 1 ) 
(A- 11) 


(39) 
(41) 

(40) 


(A-2) 


C 

12 


Vf5)  =— .J72760F-08 
IF  (TV-9000.)  45*  45,  43 
V(I)  =-.091110 
V(2)  =  .00005807 

V(3)  =  .137I74E-08 

COP  =  VI  1  I  +  TV  *  I  V(  2 )  +  TV  »  V  ( 3 1 
GL  =  COP/ ( COP-. 068 5 73? 6 )  ' 

CPREF  =  COP 

PRRFF  r  URFP  *  CPRFF  /  CRFF  *  r.A 
IF  fl-l )  50,  48,  50 
CPS  =  COP 

TV  =  J^446*#pS*SM.y  *GL**.5*  (  TS+ 198.6 )  /T 8**7 

IF  (TROOY-l.O)  50,  7o,  50 
CONTINUE 

oLA-  «L»IF  °R  C0NE  ""  LAMINAP  CONSTANTS 

CC  =  .33? 

FX  =  0,5 

IF  (RF-RET)  60.  60,  5? 

RLA=  0.892  °R  C°NE  ~-rUReULENT  CONSTANTS 
CC  =  .01396 
FX  =  0.85 


IF  ( TR0DY-3 . )  60,  54,  60 

UV  -  »  . 


HX  = 

HL-TMttIr,  ?Ror0NF  “  HEATI*r>  COEFFICIENTS 
ri  -  (1TS*UL/  R6  }**.5*PS*SM/TREF)**F* 

H  *  CC*  GA  *HX  WIURFf/y  «  ’  _ 

GO  TO  79  U.  EX)  *  v.PS/PRREF*«. 66666667 

RE  =  0.0 

I F ( VFt-6000. )  78,  76,  76 

sr?"  FT/5fc 

GO  TO  80 

SPnFRF  —  VFLOCITY  LF^S  THAN  6000  FT/SEC 
w  "  <85 

TPrr  :K:J}  *  ?2.! *  IVFL*PS/RG  ,  <*.5  /  SORO 


{  V-u  L  *  k  S , 

TR  C  =  TS* ( 1 »+.5* (GL-1 » )*R*SMSO) 

OW  =  H  *  ITDCr^rui 


OW  =  H  *  {TREC-rW!  ' 

OOOT  =  Qw-E*  I TW/ 100.  )  4 809. SF-05 

M  =  1 

IF  (N-l)  800,  9 4,  son 
LAYER  TEMPFRATURF  loop 
00  92  L=1 ,Nt 

IF  ( L- 1 )  83,  81 ,  83 
WALL  TEMPERATURE  INCREMENT 

??f  5Ti -1  f ‘SiC  ‘S^bS1  *B  cmi  * ( ^  t-reMPL  *  1 , , , 

THIN  wall  TEMPFRATURF  INCREMENT 
OTFMp  =  TEMP).  (  3  )  -  fF MPL  ( 1  ) 

OTEMP  =  {TAn<?MODOT+r{7,sr,fFMp,  _,r,_wri  „ 

TEMPI  I  )  =  TEMPL  (I  )+OTFMP  '(*^  r ( ' > /n <? I +TA  HPI/OUl 


-  to»^or*V^Ah?Sfeiv'<2fe.- 


•’awsfil 


L  =  2 
GO  TO  89 

02  =  IfMPL(L-l)  -  TFMPL  (  L ) 

END  POINT  TEMPERATURE  INCREMENT 
DTFMP  =  2,*8<M>*D? 

IF(L-NL)  84.  89.  84 
Di  =  TEMPL(L+1)-TEMPL(L> 

INTERIOR  TEMPERATURE  INCREMENT 
DTEMP  =  R ( M )  *  (01+D2) 

IF  {  NL1+1-L  I  86.  85,  89 
M  =  2 
GO  TO  88 

IE  (  NL1+NL2+1-L)  89,  87,  89 
M  =8 

INTERFACE  TEMPERATURE  INCREMENT 

DTEMP  =  2.  «  !(  D*  4  F « M  )  +D2  1  /  (  l./BIM-l  J+F(K)/B(MJ )  ) 

IF  (NL-l-L)  89,  680,  89 

INTERFACE  AND  END  POINT  TEMp  INCREMENT  Frit  THIN  INNER  MATER  1 
DTEMP  =  ( 1 • /B (M-l )  +  2,#E ( M } /B t  MI ) 

TEMP ( L )  =  TFMPL I L )  +  DTEMP 
L  =  NL 

LOCAL  TEMPERATURE 

TEMP(L)  =  TEMPL ( L ) +DTEMP 

IE(TEMP(L)-TMELT)  92.  9C0,  Q00 

TEMP ( L I  =  TMPLT 

IF.  ( L-l  )  92.  90,  9? 

MELTING  POINT  FOP  WALL 
OAR  =  OAP+ODOT  *DT I MF 
TAUAB  =  12,*0A3/(QE*RH0(1 )) 

IF  (ABOUT}  92,  91,  92 
PUNCH  1019 
PUNCH  1006,  N 

PUNCH  1005',  TIME,  ALT,  VFL 

ABOUT  =  1,0 

CONTINUE 

DO  920  L  =1 »NL 

TEMPL  ( L )  =  TFMPtL) 

IE  (SPNSE  SWITCH  3)  y5 ,  921 
I F ( SENSE  SWITCH  1)  95,93 


^  xjgi 


-c-.ym 

,  •  C'.fc-'V.-,  vjSQ 

- 


(A-12) 
(A-l  3) 


IF( TPONT-DTOUT ) 
TPRNV=0.0 
PUNCH  1006, 
PUNCH  1005.  T I 
PUNCH  1005,  OW 
IF(NCFIT)  97,9 


99  *9^*99 


N 

TIME, 

OW, 

97,97,96 


ALT, 

ODOT, 


RE 

TAUAB 


PUNCH  1005?  SM, 

PUNCH  1005,  EM, 

DO  98  L=1 »Nl 

PUNCH  1007.  TEMPL (L), 

T 1  ME -  TIME  OTIME 
TPPNT=TPRNT+DT IMF 
IE  (SENSE  SWITCH  3  5  2  , 
IF  (TIME-THI)  110,  !I0, 

ENO 


Input  F  ormat 


ft 


3.  Input  Comments 


I)  Eleven  Atmospheric  Properties  Data  Cards,  as  described 


in  Appendix  B 

2) 

TBODY = 

NCFIT  = 

Sphere 

1. 

2 

Flat  Plate 

2. 

5  or  0 

Cone 

3. 

5  or  0 

When  NCFIT  ^  0,  the  curves  needed  are 

Curve  I:  ALT,  ft  vs  Time,  sec  4  4 

Sphere  j 

Curve  II:  VEL,  ft/sec  vs  Time,  sec  }  FP 

Curve  III:  Ms  vs  M  or 

Cone 

Curve  IV:  P„/P  vs  M  I 

3  00  oo  I 

Curve  V:  T  /T  vs  M  ' 

s  oo  oo 

3}  Sphere:  Q  RN  QR  J 

FP  or  Cone:  9  Y  RET 

(9  is  used  for  identification  only  for  flat  plate  or  cone) 

4)  Melting  temperature  effective  heat  of  ablation  and 
emissivity  of  wail  material. 

5)  Number  of  materials  =  1,  2,  or  3. 

6}  Material  properties,  where  m  =  1  to  NMAT. 

Sum  of  layers  cannot  exceed  14.  NLAY  cannot  =  2. 

7)  Time  inputs:  By  card  -  SW  2  ON.  By  typewriter  on  4 
lines  -  SW  2  OFF 


NT  =  Total  number  of  time  steps 
NT  =  iT¥I-TLO)/DTIME  +  1. 


8}  Initial  temperatures  for  each  layer  and  end  point.  LI  is 
input  as  convenience  to  use  output  TEMP  cards  to  restart  -  does  not 
necessarily  equal  L. 

9)  Number  of  segments  of  the  curve.  Total  number  of  seg¬ 
ments  for  all  curves  cannot  exceed  15. 

10)  Limits  of  segment: 

X  =  Time,  for  Curves  1  and  2 


X  =  M  ,  for  Curves  3,  4,  and  5 
00 


(A -2) 
(A  -3) 


11)  Coefficients  of  the  normalized  cubic  equation: 


where 


and 


Vi  =  A0  +  A,  (FX)  4-  A2  (FX)2  +  A3'FX)3 

(A-4) 

Fv  =  {X  -  XLO)/(XHI  -  XLO) 

(A~5) 

Vj  =  ALT 

(A-6) 

V2  =  VEL 

(A-7) 

v3  =  ms 

(A -8) 

V,=  Ps/Poo 

(A-9) 

V5=  Tg/T* 

(A-10) 

12)  Constant  local  conditions,  given  for  flat  plate  or  cone  when 
NCFIT  =  0.  ALT  and  VEL  are  used  for  identification  only. 

4.  IBM  1620  Operating  Instructions 

a.  Compiling  and  Starting 

1)  Compile  Fortran  program  on  IBM  1620  40K  digital  com 
puter  starting  at  08200  memory  core. 

2)  Load  object  deck.  Check  console  switches. 

b.  Console  Switches 

1)  SW1  ON  For  output  at  all  time  steps 
*  2)  OFF  For  output  only  at  time  steps  determined  by 

DTOUT 

3)  SW2  ON  Card  input  of  time  values --D TIME,  DTOUT, 

TLO,  THI 

4)  OFF  Typewriter  input  of  these  values 

5)  SW3  ON  Ends  run,  prints  last  time  step,  branches  to 

next  case 

6)  SW4  OFF  During  typewriter  input 

7)  ON  To  make  corrections  in  typewriter  input 


-  •  r--.  -v!£^%a£&s® 


c'  Typewriter  Output  and  Input 


The  time  increment,  At  or  DTIME,  is  calculated  for  each 
material  and  pointed.  Operator  determines  At  from  smaller  value,  then 
inputs  four  time  values  by  typewriter--unless  SW2  is  ON.  The  computed 
values  for  DTIME  should  be  approximately  the  same  for  all  materials 
used.  DTOUT  must  be  an  e>act  multiple  of  DTIME. 


■v'  i 


A  rough  estimate  of  th^  machine  time  which  will  be  required  is 
printed.  This  enables  the  operator  to  leave  the  machine  and  plan  for 
additional  machine  time,  if  necessary. 


\  ~  T1620  =  NT  x  NL/40  ~  (A-ll) 

Typewriter  input  anl  output  for  a  sphere  (Example  1)  follows: 

SWi  OUTS  ALL  STEPS- -SW2  CARD  INPUTS  4  TIME  VALS 
SW3  ENDS  RUN 

TEMPERATURE  VS  TIME  FOR  LAYERS  OF  SPHERE 


DTIME  FOR  EACH  MAT 
274.  09802E-04  1 

300.54519E-04  2 

317. 39871E-Q4  3 

INPUT --DTIME,  DTOUT,  TLO,  THI--ON  4  LINES  BY  TY 
.  025RS 
.  5RS 
16.  RS 
18.  RS 

APP  1620  MIN  REQUIRED 
.  13000000E+02 


d.  After  Ablation 


K 

f 


When  the  wall  temperr.ture  reaches  the  given  melting  tem¬ 
perature,  TMELT,  ablation  begins.  After  this  point  the  computed  values 
for  the  temperatures  of  each  layer  may  be  doubtful  since  only  a  simple 
procedure  is  included  for  this  ablation. 

Qab  =  2  (Q  At),  while  Tw  £  TMELT  (A- 12) 

TAB  “  12  QAB  ^  (Qeff  Pi)  (A-13) 


30 


sowaraoee 
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Ip 


iH 
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r  -  -  *'3 
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Terminating  and  Restartins 


The  problem  may  be  terminated  before  reaching  the  upper 
time  limit,  TH1,  and  started  again  later  as  follows: 


1)  To  terminate:  Turn  SW1  ON  to  putput  values  at  all  time 
steps --to  determine  suitable  TIME  to  restart.  (SW3  on  outputs  last  time 
values  only  then  branches  to  next  case)  Save  the  cards  for  temperatures 
at  each  layer. 

2)  To  begin  again:  The  original  input  values  are  used 
except  TLO  =  TIME  and  initial  temperatures  for  each  layer  and  end  point 
are  from  last  time  step  computed.  The  time  values:  DTIME,  DTOUT, 
TLO,  and  TH1  may  be  input  on  one  card,  instead  of  by  typewriter,  if 
SW2  is  ON. 


Example  Runs  for  Sphere,  Flat  Plate.and  Cone 


Case  1  is  a  sphere  consisting  of  three  materials:  beryllium 
(4  layers),  molybdenum  {3  layers)  and  a  thin  inner  wall  of  aluminum. 

The  nose  radius  is  one  foot  and  the  position  angle  is  zero  degrees.  An 
initial  time  of  16  seconds  is  given  to  utilize  all  heat  coefficient 
equations --those  for  velocities  less  and  greater  than  6000  ft/sec.  Type¬ 
writer  input  is  used  for  DTIME,  DTOUT,  TLO  and  THI  for  Cass  1  only. 


Case  2  is  a  4°  flat  plate,  with  Y  =  2  feet,  and  composed  of  a  thin 
wall  of  beryllium  and  a  slab  of  molybdenum.  Altitude  and  velocity  data 
are  the  same  as  for  Case  1.  With  initial  time  of  16  seconds  and  transi¬ 
tion  Reynclds  number  of  22500000  both  laminar  and  turbulent  boundary 
layer  equations  are  used. 


Case  3  is  a  4°  cone  with  all  input  data  the  same  as  for  Case  2, 
except  no  curve  fit  data  are  used.  Constant  altitude,  velocity  and  local 
properties  are  given. 


Input  Data  for  Examples 

(Switch  1  on  for  first  steps) 
(Switch  2  on  for  Cases  2  and  3) 


0. 

288.  16 

101325. 

AP59 

1 

11000. 

216.  66 

22632. 

AP59 

2 

25000. 

216.66 

2488. 6 

AP59 

3 

47000. 

282.66 

120.444 

AP59 

4 

53000. 

282.  66 

58.  3215' 

AP59 

5 

79000. 

165.66 

1.  00946 

AP59 

6 

90000. 

165.  66 

. 1 04438 

AP59 

7 

31 


- ••vsrrrs* 

’*grm 

■z  i  -  -^%^’Z 


^v41gl 

105000. 

225. 66 

.  7452&5E-C2 

AP59  8 

160000. 

1325.  66 

.  362003E-03 

AP59  9 

-  -v 

170000. 

1425.  66 

.  2S2362E«03 

AP59  10 

700000. 

3325.  66 

. 000000 

AP59  11 

1.  0 

SPHERE 

0.0 

1.0 

1.  0 

&  | 

2805. 

3000. 

0.  S 

■  : 

'  'S  *■ 

3 

. 0134” 

114.  9 

.  6968 

.  032 

4 

BERYL 

i  ‘  .V,  ‘ 

. 01747 

636.8 

.  0634 

.0153 

3 

MOLY 

l 

. 02333 

168.  6 

.  2440 

.  006 

i 

ALUM 

;  _  ‘ 

540. 

1 

^  * 

540. 

& 

X 

-  /<V'V 

540 

1 

:  -  Y 

.  *1  ' 

540. 

1 

^ j  „  *.rC 

v  .  r  ^ 

.5  » -x 

540. 

1 

y  1* 

540. 

1 

540. 

1 

. 

540. 

1 

540. 

1 

1 

.. 

0.0 

• 

00 

rH 

300. 

1 

0.0 

-255.25706 

48. 9861 

-93.729 

18. 

V  ' 

16000. 

4370. 0192 

-16559.985 

-810.016 

i  ' 

2.  0 

5 

FLAT  P 

J 

4.0 

2. 

22500000. 

2805. 

3000. 

0.  8 

2 

' 

.01347 

114.9 

.  6968 

.  008 

1 

BERYL 

. 01747 

636.  8 

.  0634 

.  0153 

3 

MOLY 

.  025 

.  50 

16. 

18. 

540. 

1 

540. 

1 

> 

540. 

1 

540. 

1 

540. 

1 

0.  0 

1 

18. 

300. 

* 

-255. 25706 

48. 9861 

-93.729 

t 

0.0 

18. 

16000. 

4370. 0192 

-16559. 985 

-810. 016 

32 


1.  5 

1.4690000 
3.  0 
2.92 
6.  0 
5.689 
10. 

9.065 

4 

1.5 
1.  045 
3.  0 
1.  13 
6.  0 
1.  39 
10. 

1.947 

4 

1.5 

1, 01231 
3.0 

1. 0351800 

6.  0 

1.0987000 

10. 

1. 2092 
3.0 

4.0 

2805. 

2 

. 01347 
.  01747 
.  025 

540. 

540. 

540. 

540. 

540. 

45.  98 
6. 076376 


1.4689996 

-. 013494600 

-. 00450200 

6.0 

2, 8634991 

-. 09449200 

-. 00000500 

10. 

3, 5613287 

-.  18000000 

-. 00532700 

20. 

7. 9674930 

-1. 2445000 

. 03500000 

3.  0 

. 07599960 

-. 02249760 

, 03149860 

6.  0 

. 19699980 

. 09450330 

-. 03150230 

10. 

.49766601 

-. 05000000 

.  10933290 

20. 

1. 5800819 

. 95175000 

-.  13583300 

3.  0 

. 01994460 

. 00153250 

. 00139350 

6.  0 

. 05312460 

. 01116250 

00076650 

10. 

. 09385967 

. 0208000 

00416070 

20. 

. 31294110 

. 09418250 

. 00802500 

0 

2. 

2250C000. 

3000. 

0.  8 

114.  9 

.  6968 

.  008 

636.  8 

.  0634 

.  0153 

.  50 

1 6. 

18, 

1 

1 

1 

1 

5 

6231. 
425,  852 


432. 5785 


CONE 


BERYL 

MOLY 


Output  Data  for  Examples 


a.  Temperature  versus  Time  for  Layers  of  Sphere 


i hFTA=  RN-  QR= 

,OOO000O0F-99  .lQOOOOOOE+OI  .300O0000E+O1 

tmfi.t=  OEFF=  EMISS1VITY= 

• 28050000R+04  .30000000E+O4  .80000000E-00 

MATERIAL  —  OTIME.  H»  K *  RH0.  CP»  TAUT,  NlAY 
274.09802F-04  3 

. 3 3470000F-03  .  1 IA90000F+O3  .69680000F-00  .80OO0O0OF-02 

30O.54519F-04  2 

. 3 7470000F-01  .63680000E+03  .63400000E-01  ,51000000F-02 

337 . 39871 F-04  3 

.23330OO0E-01  • I6860000F+03  .24A00000F-00  .60O00000E-02 

DT IME=  OTOUT  *  TLO=  THI  = 

•25000000F-01  . 50000000F-00  . 3 6000000F+O2  , 1 7500000F+0? 


**-  FOR  FACH  OUTPUT  TIMF  STEP  - 

TIME  ALT  »KFT  VEL,  FT/SEC 

QW  QDOT  H  TAUAR  - 

MACH  NO.  PRESSURE  TEMPFRATURE  - 
MACH  NO.  PRESSUPF  TFMPFRATUPF  - 


RF 

-  HF  AT  COERES 

-  LOCAL 

-  FRFF  STREAM 


TEMPERATURE  AND  COUNT  FOR  EACH  SUN  LAYER 


. 36000000F+02 
• 83836260F+02 
.000000O0F-99 
•6A375I82E+01 
5AO.000O0F-0O 
540, 00000 F-0O 
5AO,00OO0P-00 
5Ae.n0OO/^F-0'' 
5AO,r»0On0F-00 
s^'^coonoc-ne 
5AO.O0000F-0O 
540,oOooof-oo 
54O.00OO0F-0O 
2 

*  16O25OOOF+02 
.83*33656F+0? 
,  0enrH;.''00F-9Q 
.6398  753AF+0  3 
546.505  3  3  F-00 

540,O0O00F-0n 
5AO,n0O00F-0O 
54O.O0O00F-OO 
5A0 . OOOOOr - 00 
540  ,e00O0r_0f> 
54O.O0O00F-0O 

5i,n,noooor-0'' 

5/.rt,fi0O0OF-0O 

A 


•A598I 030F+O2 
,83832°89P+02 
.  3  587273  AE  +  05 
.29A91558E+03 
3 
2 

3 

4 

5 

6 

7 

8 
*9 

. A5A385 l OF  +02 
•  83330385F  +02 
.36097295F+05 
.302693  A7F  +  0  3 
1 
2 
3 
A 


.6231 1 290E+04  .00000000^-99 

.OOOOOOOOF-99  . 00000000^-99 
.3622 3369^+04 
•  38998800E+O3 


•61936070F+OA  . OOOOOOOOF-99 

.0000000^^-99  ,i'0OO00O0c-99 

.3SP36250F+04 
. 389988 00F  + 0  3 


34 


7 

. 16150000F+02  .42714160F+02 

• 79992499F+02  . 79988578E+02 

• OOOOOOOOE-99  . 17253935E+05 
.62038092F+01  .34495806F+03 

568,40264F~00  1 

544.48412F-00  2 

54O.41221F-00  3 

540.02213F-00  4 

54O.00094E-CO  5 

540.00O07F-0O  6 

540.00000F-00  7 

540.00000F-00  8 

540.00000E-00  9 

8 

.16175000P+02  .42166950E+02 

• 88532884F+0?  •65528566F+G2 

• OOOOOOOOF-99  . 1 749 183 7F+05 
•61645947F+01  .35413596E+03 

572* 16971F-00  1 

545.78844F-00  2 

540 .654 18 F— 00  3 

540.04637F-00  4 

54O.O0287F-0O  5 

540.00040F— 00  6 

54O.00OO2E-00  7 

54O.00OOOF-0O  8 

54O.00000F-0O  9 

9 

• 16200000F+02  .41618950E+02 

. 87703263E+02  .87699I40E+02 

.OOOOOOOOF-99  • 177319P6F+05 
.61253071F+01  . 36357234E+03 

575.54830F-0O  1 

547. 1 8479 F— 00  2 

54O.95165F-00  3 

540.08345F-00  4 

540.00669F-00  5 

54O.O0126F-0O  6 

540.00016F-00  7 

54O.00O00F-0n  8 

540 . 00O002  90  9 

10 

. 16225000F+02  .41070170E+02 

•86877063F+02  .868 7284 1 F+02 

.OOOOOOOOF-99  .17974121F+05 
•60859443E+01  . 37327468F+03 

578.60184F-0O  1 

548.63920F-0''  2 

541. 30423F-00  3 

540.13546F-00  4 

54O.01322F-0O  5 

540.00306E-00  6 

540. 0005 5F -00  7 

54O.00O03F-00  8 


.600491 30r +04  .OOOOOOOOF-99 
e OOOOOOOOF-99  .OOOOOOOOF-99 
.339 19010E+04 
.38998800F+03 


.F9669560F+04  .OOOOOOOOF-99 
•31781561F-01  s OOOOOOOOF-99 
.33540708F+04 
•38998800E+03 


.59289230F+04  .00900000E-95 

•31959027E-01  .OOOOOOOOF-99 
.33164104F+04 
.38998800F+03 


•58908270E+04  .OOOOOOOOE-99 
.321 36568E-01  .OOOOOOOOE-99 
. 32  7891 92F  +  04 
.38998800E+0^ 


b.  Temperature  versus  Time  for  Layers  of  Flat  Plate 


THETA=  Y  =  RET* 

.40000000E+01  « 20000000E  +01  .22500000E+08 


k  r %m 

h. 


TMELT=  0FFF  =  EMTSSIVTTY= 

♦28050000F+04  .  3Q000000E+C4  .  80000000F-RO 


MATERIAL  —  DTIMEt  M,  K»  RHCt*  CP.  TAUT »  NLAY 
274.09802F-04  1 

. 13470000E-01  • 1 14900C0E+03  .69680000E-00  .80000000E-0/ 

300.545 J9F-04  2 

. J7470000F-01  .63680000E+03  .63400000E-01  .51000000E-02 


DTIME=  DT  OUT  =  T  L0=  THI  = 

*  25000000E-01  « 30000000E-00  .16000000E+02  . 18000000E+Q2 


**-  FOR  EACH  OUTPUT  TIME  STEP  - 

TIME  ALT.KFT  VEL. FT/SEC  RE 

GW  QPOT  H  TAUAR  -  HEAT  COEFFS 

MACH  NO.  PRESSURE  TEMPFRATUPP  —  LOCAL 
MACH  NO.  PRESSURE  TEMPERATURE  —  FREE  STREAM 
TEMPERATURE  AND  COUNT  FOR  EACH  SKIN  LAYER 


1 

. 16000000E+02  .4598 1030E+02 

. 559785 12E+01  .55945795E+01 

•  60  763760E+0 1  .42585199E+03 

•  64375 182E+0 1  . 2949 1 558F+03 

54O.00000f-*00  1 

540.00000E-0D  2 

540.00000F— Or*  3 

540.00000F-00  4 

540. 00000 E-00  5 

2 

.  16025000F+02  .454385  1  OF  ■‘■02 

.55888277F+01  .55855560E+01 

.60422289F+01  .43564O50E+03 

.63987534E+01  . 30269 147C +03 

540.18782R-C0  1 

540. 18782F-00  ? 

340,0 00 OOF— 00  3 

540.00000F-00  4 

540.O0O00F-0O  5 

3 

. 16O50OO0F+02  .44895100F+02 

. 52755893F+02  .5275261 7F+02 

.60079820F+01  .44566301F+03 

.63599142F.+01  . 31068475F+03 

541 .94009E-00  1 

541 .94009F-00  2 

540,0 1 234F-0O  3 

C.40,OOOOOF-O0  4 

540,n0ooor-on  4 


•6231 1290E+04  .21965317E+08 

•22409938E-02  .00000000E-99 

.43257847E+03 
• 38998800F+03 


•61936070E+O4  .22369129E+08 

. 22638823F-02  . OOOOOOOOF-99 

.43220779F+03 
.389988G0E+03 


. 6 1 5601 30F+O4  . 22779691 F+0 8 

•20566093E-01  . OOOOOOOOE-99 

.431 83789E+03 
. 38998800F+03 


36 


m 

I 

1 

ft 

1 

4 

i 

. 160750008+0? 

.443511108+02 

.611834608+04 

•  2  3 197024  8+03 1 

£§ 

.531017026+0? 

.430983846+02 

.209671068-01 

.000000008-99 

H' 

.  597363398+0! 

.444°  ?4  *)9F  +03 

.43  14686i!8+03 

.632099928+0  l 

.318901 16F+03 

. >89988008+03 

§ 

543.5O3=7=-0o 

1 

§E 

5A3.59357P-00 

2 

44o,  )  38?  2  ='-00 

445,0908  1  F-(in 

4 

P 

840.000008-0^ 

5 

| 

. 161 00000F+0? 

.438062508+02 

.60806G7O8+04 

.236212208+0.4 

R 

. 534501 318+02 

« 5  34467748  +02 

.’13770288-01 

.000000006-99 

lit 

. 5939 1 8608+0 1 

.466429788+03 

.431 100226+03 

& 

.623201006+01 

.3  >7147546+03 

.389588008+03 

l 

644.9131 '  =  -00 

1 

644.91 31 7p— 00 

2 

1 

640,35n?6F— 00 

? 

g- 

5  4  0  *  0  9  7  0  f.  ■**  0  0 

4 

540 .000 1  OF— OO 

s 

'$ 

'  1 

. 16500000K+0 2 

.349811008+0? 

.546694508+04 

.303780298+08 

if-- 

. 885?069 1 F+0? 

.  485 16958K+02 

.•>9  1020458-01 

,OOoOOOOO=-99 

.  534  70380F  +  0 1 

.673814016+03 

.429438368+03 

p 

. 66181 3476+0 1 

.4980291 16+03 

. 394  I  43448  +  03 

6  48 . 7467  36-00 

1 

1- 

=38.746738-00 

2 

64  7 .05337F-00 

H 

842.?0022c-00 

4 

l 

441 .00586F-OO 

5 

' 

n 

♦ i7000COOF+n? 

.236589908+02 

.467376806+04 

.35624770?  +o« 

J& 

.597666476+02 

.597626198+0 2 

. 4045 1 9=?F-01 

.000.300008-99 

.439731776+01 

.104337266+04 

.463098808+03 

H 

.457502798+01 

.831427946+03 

.434411698+03 

||- 

569 . 2  7844  =— 90 

1 

Wr 

669.278546-00 

2 

P 

6  66 ,628°.  9F -00 

If 

846. 1  34668-00 

4 

H 

f 45.30388F-0O 

5 

s 

61 

p 

. 175000008+02 

.120026706+02 

.385149408+04 

.389275738+0“ 

8 

. 53IQ2654F+02 

.5? 18841 7?+0? 

.529597C5F-0! 

.000000006-99 

p 

. 549081046+01 

. 15773304= +04 

.497934926+03 

8 

.36O20070F+01 

. 134460’0F+O4 

.475907O9C+03 

jl 

4  76  .  ''23  =  48-00 

1 

ft 

576  » 3 234 4P— 00 

2 

H 

563. 16318r-Oo 

3 

at 

545. 1 4646F-00 

4 

p 

4* 2, 483278-00 

fc 

81 

. 1RO00OO0F+02 

.400000006-04 

.30000] 806+’4 

. 3967987~6  +08 

& 

. 3006  7380F  +  0  ? 

.390430426+02 

.441  10878F-O1 

3  nooO0000c-QO 

8 

•  26?  1-26  =  98  +  0  1 

.234199028+04 

.434  •  1 9408  +  03 

|g 

.26874580? +Q1 

.21 1 67 1 058  +04 

.41  S687P6"  +  >'3 

M 

47'}.  5b?8 1  F-0O 

1 

n 

=  79.4628 1F-0O 

2 

s 

568.944410-00 

3 

1 

461.828518-00 

4 

11 

=45.349208-00 

c, 

TH6T  A  =  Y  =  RET  = 

•  0000006+01  .200000006+01  .  225000006+08 

TMFLT=  OEcF  =  EMI  SSIVITY® 

.289500008+04  .300000006+04  .800000006-00 

MATERIAL  —  0TIM6,  M,  K.  RH0»  CP.  TAUT,  NLAY 
274.09802F-04  1 

.134700006-01  .114900006+03  .69680000E-00  .80000000F-02  1 

300.945 196-04  2 

.174700006-01  .63680000E  +0  3  .604000006-01  .51OOO00OE-02  3 

DTIME=  OTOUT  =  T  L0=  THI  = 

.25OO00O0F-01  ,800000006-00  .160000006+02  .180000006+02 

MACH  NO.  PRESSURE  TEMPFRATURF  —  LOCAL 
.60763760F+01  .425852006 +03  .43257850^+03 

**-  FOR  EACH  OUTPUT  TIMF  STEP  - 

1  IMF  ALT,<?t  VEL.FT/S6C  RE 

Ow  ODOT  H  TAUAR  -  HEAT  COFFFS 

1+MPERATURE  AND  COUNT  FOR  EACH  SKIN  LAYER 

1 

. 16O00000F+02  .459800006+02  .62310000F+04  .219653156+08 

.96957630F+01  .9692  71 3F+01  .308151496-02  .909000006-99 

540.00000F-00  1 

540.000006-00  2 

540.00000^-00  3 

64O.00OOOF-0O  4 

•^O.OOOOOF— 00  5 

2 

. 160250006+02  .459800006+02  .623100006+04  .219653156+08 

.969576306+01  .969249136+01  .388151496-02  .000000006-99 

540.325926- 09  1 

540.325926- 00  2 

540.00O00F-0O  3 

540.090006— Oo  4 

540.00O09C-00  5 

3 

, 16O500O0F+05  .4598000 OF +02  .62310000^+04  . 2 1 9653 1 56+03 

. 96939 1 4 i 6+0 l  .9690634  56+01  ,388 1 36586-02  . OOPOOOOOC-99 

540.614746-00  1 

‘40.61*246-00  2 

54O..9?!4lr-00  3 


• 

81 

.  180000.90F+02  .4  5980000F  +0  2  .6’3  1 0O00F  +  04  .  2  1  965  3  1  56+08 

.-'-65353146  +  01  .  6504/736  +01  .38  7012586-''/  .OOOOOOOOE-99 

,  oioononoc-ci?  .000000006-99  .ono^ponop-og  .74000000F-00 
.  i  I  !2b?70F  +  04  .(•-4262-’«6>r-06  .  +  ?0,n54656-O6  .  24  73  1  OO-.F-00 

+47.436838—09  ] 

847, 4'f 8TT-0O  ? 

>  5  .  ]  4  7820 -00  » 

84’.70!,'3r-9n 
5-i  3. 325766-00 


4 

5 


Appendix  B 

1559  ATMOSPHERIC  PROPERTIES 

Table  I  presents  1959  atmospheric  properties  used  as  input  to  the 
Fortran  program. 


Table  I.  ARDC  1959  Atmospheric  Properties 


ALT,  Meters 

■m 

P,  Newtons /m2 

Card  # 

0.  0 

288. 16 

-  - 

101325. 

1 

11000. 

216.  66 

22632. 

2 

25000. 

216.  66 

2488.  6 

3 

47000. 

282. 66 

120.444 

4 

53000. 

232. 66 

58. 3215 

5 

79000. 

165. 66 

1.  00946 

6 

90000. 

165. 66 

.  104438 

7 

105000. 

225.66 

. 00745265 

3 

160000. 

1325. 66 

.362003  E-03 

Q 

t 

170000. 

1425. 66 

.282362  E-03 

10 

700000. 

3325. 66 

0.0 

11 

The  following  changes  in  dimensions  are  made,  and  the  properties 
are  stored  into  memory  as: 


AP  (K,  1)  =  ALTft  =  ALTmeters  X  3.  2808333 
AP  (K,  2)  =  T,  °R  =  T,  °K  X  1.8 


(B-t) 

(B-2) 


AP  (K,  3}  =  P,  lbs /ft2  =  P,  Newtons/m2  X  .0208854  (B-3) 


where  K  =  1  to  11 


The  local  velocity,  V<jo,  is  given  in  ft/ sec  and  the  local  altitude, 
ALT,  i  a  given  in  kilo-feet  (geo  metric  measure)  and  is  converted  to  feet 
(geopotential  measure). 

A.LTF  =  20856000.  (ALT)  /  (20856.  +  4LT)  (B-4) 

AALT  =  DALT  *  ALTF  -  APK  (B-5) 

»  l 

Too  =  APk>  2  +  ALT  (  TS)  (B-6) 


where  TS  is  the  temperature  slope  for  the  atmosphere  layer. 


When  TS  -  0.  0 


where 


and 


When  TS  £  0.  0 


where 


Then 


and 


_  A  AP  Temperature 
A  AP  Altitude 


Poo  =  APk,  3/  eJ 


e  =  2.  718281828 


j  =  .  01879  AALT/APK|  2 


Poo  =  APk>  j 


APk, 


L  T‘o 


j  =  .  01879/TS 


Poo  =  .01879  Poo/T^ 


^oo  ~  Voo  l^/Rg  T,  o 


(B-7) 

(B-8) 

(B-9) 

(B-10) 

(B-l  1) 

(B- 12) 

(B-13) 

(B- 14) 
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SYMBOLS 


TBODY 

NCFIT 

TH 

RN 

Y 

RNORY 

QR 

RET 

RE 

QRRET 

TMELT 

QE 

E 

NMAT 

C(M) 

RHO(M) 

CP(M) 


0 

D 

R 

Y 


c  /  q 

Ret 

Re 


Teff 


Pm 


C 


Pm 


Type  of  configuration  -  1.  Sphere 

2.  Flat  plate  3.  none 

Number  of  curves  fitted  for  input  data 

Angle,  deg:  Sphere;  position  angle 
Flat  plate  or  cone;  used  for  identification 

Nose  diameter  of  sphere,  ft 

Nose  radius  of  sphere,  ft 

Length  along  surface,  ft 

Input  symbol  (RN  or  Y) 

Ratio  of  laminar  heat  transfer  to  stagnation 
rate  heat  transfer,  from  Figure  5 

Transition  Reynolds  number,  given  for  FP 
or  cone 


Local  Reynolds  number,  defined  by  Equation 
(36)  for  FP  or  cone 

Input  symbol  (QR  or  RET) 


Melting  temperature  for  wall,  °R 
Effective  heat  of  ablation,  Btu/lb 


Emissivity  of  wall  material 
Total  number  of  materials;  1,  2  or  3 
Material  thermal  conductivity,  Btu/ft-sec  -  °R 
Material  density,  lbs/ft 

Specific  heat  for  the  material,  Btu/lb  -  °R 


TAUT 

Ttot 

Total  thickness  of  the  material,  ft 

TAU(M) 

Tm 

Thickness  of  each  layer  of  the  material,  ft 

NL1  -2-2 

NLA.Ym 

Number  of  layer  s  of  each  material 
(Total  cannot  exceed  14) 

F»  3 

F  ,  B 
rn'  m 

Functions  of  material  properties,  defined  by 
Equations  (39)  to  (4i) 

TIME 

t 

.Local  time,  sec' 

DTIME 

At 

Time  increment,  sec,  used  for  calculations, 
Equatior  (38) 

DTOUT 

Time  in:remen':,  sec,  used  for  output; 
Multiple  of  DTIME 

TLO 

41 0 

Initial  time,  sec,  TLO  ^  0.  0 

THI 

^i 

Upper  time  limit,  sec 

NT 

Total  number  of  time  steps,  Equation  (A-l) 

T1620 

Estimate  of  1620  machine  time  required, 
Equation  (A-l  1) 

t 

TPRNT 

Count  for  time  print-out 

ABOUT 

Test  for  output  of  TIME,  ALT,  and  VEL  at 
first  ablation 

NUMj 

NSEG 

Number  of  segments  to  each  of  2  or  5  curves 
given  (Total  number  of  segments  cannot 
exceed  15) 

CONlt  2 

Xlo, 

Xhi 

Limits  of  curve  segment,  input  data 

con3> 

4,  5,6’ 

Ao, 

1,2,3 

Normalized  cubic  curve  fit  constants,  used  in 
Equation  (A-4) 

T  Temperature  for  each  skin  layer  at  each  time 

step,  ’R;  given  for  initial  time  step,  then 
defined  by  Equations  (49)  to  (55) 


TEMP 


TEMPL 

T' 

Temperature  for  each  layer  at  previous 
time  steps 

DTEMP 

AT 

Temperature  increment  for  the  local  layer, 
Equation  (42) 

TW 

T 

X. 

Wall  temperature,  °R,  defined  by  Equations 
(49)  and  (53) 

TEMP(NL) 

tnl 

Temperature  for  innermost  point,  defined  by 
Equations  (52)  and  (55) 

D1  -Z-j 

Dl,  2,  3 

Incremental  temperature  distribution,  defined 
by  Equations  (43)  to  (45) 

FM 

Moo 

Free  stream  Mach  number.  Equation  (B-14) 

SM 

Ms 

Local  Mach  number 

Sphere:  from  Equation  (25) 

FP  or  Cone:  given  constant  or  f(Moo), 

Equation  (A -8) 

PF 

Poo 

Free  stream  pressure,  lbs/ftz,  Equations 
(B-8)  and  (B-ll) 

PS 

Ps 

Local  pressure,  lbs/ftz 

Sphere:  from  Equation  (23) 

FP  or  Cone:  given  constant  or  f(Moo), 

Equation  (A -9) 

PSYPT 

Ps/Ptot 

Ratio  of  local  pressure  to  stagnation  pressure. 
Equation  (24) 

TF 

Too 

Free  stream  temperature,  °R,  Equation  (B-6) 

TS 

T 

A  s 

Local  temperature,  PR 

Sphere:  from  Equation  (27) 

FP  or  Cone:  given  constant  or  f  (Meo), 
Equation  (A-10) 


AP 


APX,  i 


Atmospheric  Properties  of  Table  I  in  Appendix 

Altitude,  ft,  geopotential  measure, 

Equation  (B - 1 ) 


Atmospheric  temperature,  °R,  Equation  (B-2) 
Atmospheric  pressure,  lbs/ft2,  Equation  (B -3) 


APK,  2 

Atmospheric  temperature,  R,  Equation  (B-2) 

APK,  3 

Atmospheric  pressure,  lbs/ft2,  Equation  (B -3) 

ALT 

ALT 

Local  altitude,  kilo-feet,  geometric  measure, 
Equation  (A -6) 

AI-TF 

Local  altitude,  ft,  geopotential  measure, 
Equation  (B  -4) 

DAI  T 

A  ALT 

Difference  between  local  and  layer  base  altitude, 
ft,  geopotential  measure.  Equation  (B-5) 

VEL 

Local  free  stream  velocity,  ft/sec 

Equation  (A -7) 

DF 

Poo 

Free  stream  density,  ibs/ft'.  Equation  (B-13) 

p  sea 

level 

Sea  level  density,  0.  002378  -slugs*' ft3 

QAB 

qab 

Heating  rate  during  ablation.  Equation  (A-12) 

TAUAB 

Tab 

Preliminary  estimate  of  total  ablation  thickness, 
in.  ,  Equation  (A- 13) 

QW 

qw 

Aerodynamic  heating  rate,  Btu/ft2  -  sec 
Equations  (28)  and  (47) 

QDOT 

4 

Defined  by  Equation  (48; 

TREC 

Recovery  temp ~ii  ature,  °R,  Equation  (46) 

H 

H 

Aerodynamic  heat  coefficient,  Btu/ft2  -  sec  -  "R, 
defined  by  Equations  (15),  (19),  (29),  (32), 
and  (35J 

Z 1-2-3 

Z  i  j  3 

Factors  o*  Equations  (7),  (15)  and  (19).  Defined 
by  Equations  (8)  to  (10)  and  (20)  to  (22) 

HX 

Flat  plate  to  cone  heat  transfer  factor,  as  used 

in  Equations  (32)  and  (33) 

h  Enthalpy,  ft2 /sec2,  used  in  Equation  (28) 
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CHEF 


CREF 

k* 

Reference  thermal  conductivity,  Btu/ft-sec-°R 
Equations  ( 1 2a)  and  ( 1 2b) 

UREF 

M-* 

Reference  viscosity,  lbs-sec/ft2,  Equation  (13) 

PRREF 

Pr* 

Reference  P'ranatl  numher,  Equation  (14) 

TREF 

T* 

Reference  temperatui  e,  ®R,  Equation  (11) 

COP 

CP 

Specific  heat  for  air.  Equation  (56) 

G 

V 

Ratio  of  specific  heats;  for  air,  y  =  1.4 

GL 

Ys 

Computed  Y  as  a  function  of  temperature. 
Equation  (57) 

R 

R 

Recovery  factor:  laminar,  R  =  0.  85,  turbulent, 
R  =  0.892 

GA 

ga 

Acceleration  for  gravity;  32.  174  ft/secz 

RG 

Rg 

Gas  constant,  for  air,  Rg  =  1716  ft2  /  sec2  -  °R 

<T 

Stefan-Bollzmann  constant,  0.48096  X  10“5 

Btu 

ft2-sec-°R4 

St 

Stanton  number.  Equations  (1)  and  (17) 

Nu 

Nusselt  number 

Cf 

Skin  friction  coefficient,  Equations  (3)  and  (18) 

Leading  edge  Cf  proportionality  factor  for 
laminar  and  turbulent  boundaries 

8 

Velocity  gradient  parameter 
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*x*em  r**£*+*^ 


AB 

co 

S 

M 

NL 

L 

N 

I 

inc 

o 

TOT 

K 

FP 

VT 


SUBSCRIPTS 

Ablation  properties 
Free  stream  properties 
Local  properties 

Material  properties,  counts  materials  from 
1  to  NMAT 

Total  number  of  skin  layers  plus  end  point, 
NL  g  15 

Number  of  the  layer,  from  1  to  NL 
Counts  time  steps 

Counts  curves  fitted  for  input  data,  from 
1  to  NCFIT 

Incompressible  conditions 
Stagnation  conditions 
Total  conditions 

Number  of  base  of  altitude  layer 

Flat  plate 

Wall 

SUPERSCRIPTS 

Reference  properties 
Temperatures  at  previous  time  step 
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